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Abstract. We present molecular dynamics simulations on the structural relaxation 
of a simple bead-spring model for polymer blends. The introduction of a different 
monomer size induces a large time scale separation for the dynamics of the two 
components. Simulation results for a large set of observables probing density 
correlations, Rouse modes, and orientations of bond and chain end-to-end vectors, arc 
analyzed within the framework of the Mode Coupling Theory (MCT). An unusually 
large value of the exponent parameter is obtained. This feature suggests the possibility 
of an underlying higher-order MCT scenario for dynamic arrest. 

PACS numbers: 64.70.Pf, 83.80.Tc, 83.10.Rs 
1. Introduction 

Polymer blends are soft-matter systems which exhibit 'dynamic asymmetry' in the 
meaning that, starting from two homopolymers with different mobilities, two separated 
segmental dynamics can still be observed in the blend. Phenomenological approaches 
usually consider thermally driven concentration fluctuations [lj and self-concentration 
effects induced by chain connectivity [2] as key ingredients for structural relaxation in 
polymer blends [3]. A recent approach combines self-concentration effects with ideas 
of the Adam-Gibbs theory [4]. For most of the investigated systems, dynamics of the 
two components in the blend display qualitatively similar features. However, recent 
experimental results by nuclear magnetic resonance (NMR) [HIE]) dielectric spectroscopy 
[3 El El EH ED] ) or neutron scattering [121 H3] suggest that a rather different scenario 
arises when the two homopolymers exhibit very different glass transition temperatures. 
Hence, for dilute concentrations of the fast component, the two components in the 
blend exhibit strong dynamic inmiscibility. A large separation in their relaxation times 
is observed, which can be of even 12 orders of magnitude in blends of poly(ethylene 
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oxide) /poly (methyl methacrylate) (PEO/PMMA) for extreme dilution of PEO [6J. In 
such conditions the motion of the chains of the fast component takes place in a slowly 
relaxing matrix formed by the slow component, providing a connection with the problem 
of confinement in host media with interconnected voids. 

We have recently performed an investigation on the structural relaxation dynamics 
of a simple bead-spring model for polymer blends [H]. The introduction of monomer 
size disparity between the two components induces a large time scale separation for 
low concentrations of the fast component, which displays unusual relaxation features. 
Hence, density- density correlators exhibit logarithmic decays over time intervals of even 
four decades and a concave-to-convex crossover by varying the thermodynamic state 
point (i.e., the control parameters) or the wave vector [14] . Dynamic features observed 
for this simplified model are supported by recent fully atomistic simulations on the blend 
PEO/PMMA p2]. 

We have discussed the unusual features reported in Ref. [H] within the framework 
of the Mode Coupling Theory (MCT) of the glass transition [151 [IB], and suggested 
an underlying higher-order MCT transition as the origin of the observed anomalous 
relaxation scenario. Higher-order MCT transitions were initially predicted by schematic 
models [17] . and later derived for simplified models of short-ranged attractive colloids 
[T8| [15] . These systems show two different mechanisms for dynamic arrest: steric 
repulsion characteristic of colloidal systems, and formation of reversible bonds, induced 
by the short-ranged attraction. Coexistence of both mechanisms of very different 
localization lengths [T8l [T9] yields a higher-order MCT transition in a certain region 
of the temperature-density plane. The mentioned anomalous relaxation features are 
derived from the MCT equations as specific solutions associated to the higher-order 

point [H1CL91 [201 [21]. 

Results for the mean squared displacements and density-density correlators in the 
bead-spring polymer blend of Ref. [2] display striking similarities with qualitative 
features associated to higher-order MCT transitions. Similar results have also been 
observed in later simulations of binary mixtures of non-bonded particles with large 
size disparity, both for soft [221 [23] and ultrasoft interactions [23]. Finally, very 
recent two-dimensional NMR experiments on a polymer-plasticizer system have revealed 
logarithmic relaxation for the strongly confined plasticizer [25]. Hence, this collection 
of similar experimental and simulation results suggest a common relaxation scenario 
for multicomponent systems exhibiting strong dynamic asymmetry. Moreover, the 
mentioned analogies with short-ranged attractive colloids suggest that the higher-order 
MCT scenario might be a general feature of systems showing several mechanisms for 
dynamic arrest. For the mentioned polymeric and non-polymeric mixtures, we have 
suggested bulk-like caging and confinement [131 [221 [23] . These mechanisms would be 
respectively induced by the presence of neighbouring small particles and by the slow 
matrix formed by the large particles. 

It is worth mentioning that solutions of the MCT equations for a fluid of hard 
spheres confined in a disordered matrix of strictly static obstacles explicitly reveal the 
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existence of a higher-order transition |26j. As discussed in Ref. [26j, the strictly static 
nature of the matrix induces differences with the former mixtures, where the matrix 
shows a slow relaxation. Hence, though they share common features for the dynamics 
of the confined component, a comparison of results between both kind of mixtures must 
be taken with care. 

The test of MCT predictions for the bead-spring blend model reported in Ref. [H] 
was restricted to density- density correlators. In this article we present a systematic 
test for a large set of correlators probing different dynamic features as Rouse modes or 
orientations of bond and chain end-to-end vectors. Consistently with MCT predictions, 
a common set of dynamic exponents provides a good description of dynamic correlators 
in the early-middle stage of the structural a-relaxation. Consistently with previous 
results [H], the unusually large value obtained for the exponent parameter suggests 
that the observed anomalous relaxation features might be associated to an underlying 
higher-order order MCT scenario. 

The article is organized as follows. In Section 2 we summarize the main details of 
the simulated model. In Section 3 we present simulation results for static correlations. 
The main predictions of MCT are exposed in Section 4. We discuss within the framework 
of MCT relaxation features of the slow and fast component in, respectively, sections 5 
and 6. Conclusions are given in Section 7. 

2. Model and simulation details 

The model introduces a binary mixture of bead-spring chains (of the species A and 
B). Each chain consists of iV = 10 monomers of mass m — 1. All the monomers in 
a same chain belong to the same species (i.e., all them are A-like or B-like). Two 
given monomers (placed at a same chain or at different ones) interact through a soft- 
sphere potential plus a quadratic term, V Q/ g(r) = 4e[(cr Q/ g/r) 12 — C + C 2 (r/cr a/3 ) 2 ], where 
e = 1 and a, [3 G {A, B}. The interaction is zero beyond a cutoff distance r c = ca a/ 3, 
with c = 1.15. Continuity of potential and forces at r = r c is guaranteed by setting 
the values Co = 7c -12 and C*2 = 6c~ 14 . The radii of the different pair interactions 
are <tbb = 1, °"aa = 1-6<tbb, and <tab = 1-3o"bb- Chain connectivity is introduced 
by a FENE bonding potential [27], V^ NE (r) = -kR 2 eln[l - (R a aa )~ 2 r 2 }, between 
consecutive monomers, with k = 15 and Rq = 1.5. The superposition of V a p(r) and 
V^f NE (r) provides an effective bonding potential for connected monomers with a sharp 
minimum at r = 0.985a aj 3, which makes bond crossing impossible. 

The blend composition is defined as xb = N B /(N\ + N B ), where N a is the number 
of a-chains. All the data presented here correspond to a fixed composition xb = 0.3 
(we have simulated a mixture of Na = 210 and Nb = 90 chains). We use a packing 
fraction <fi = (tt/6)L^ 3 [Na0' A a + ^b°bb] = 0.53, with L the side of the cubic simulation 
cell. The value <ft = 0.53 is comparable to those used in simulations of slow relaxation 
in simple liquids [28| [29]. In the following temperature T, distance, wave vector q, and 
time t will be given, respectively, in units of e/ks, o~bb, o"bb> an d c r BB("i/e) 1 ^ 2 - 
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The system is prepared by placing the chains randomly in the simulation cell, with 
a constraint that avoids monomer core overlapping. The Newton equations of motion 
are integrated in the velocity Verlet scheme [30], with a time step ranging from 2 x 10~ 4 
to 5 x 1CT 3 for, respectively, the highest and the lowest investigated T. Standard 
periodic boundary conditions are used for calculation of monomer-monomer distances 
entering in the interactions. Computational expense for the latter calculation is reduced 
by implementing a standard link-cell method [30]. The system is thermalized at the 
selected temperature by periodic velocity rescaling. Then the equilibrium run for data 
acquisition is performed in the microcanonical ensemble (constant energy). During this 
run no drift in thermodynamic quantities is observed, either aging effects in dynamic 
correlators computed for different time origins. Statistical averages at a given state 
point are performed over typically 20-40 independent runs. 

3. Static properties 

In this section we provide information about static correlations in the bead- 
spring blend. We compute normalized partial static structure factors, S a p{q) = 
{p a (q, 0)p^(<?, Q)) / (N y/N^Np) . The quantity p a (q,t) is the density fluctuation for wave 
vector q and is defined as p a (q,t) = ^ exp[zq • r"(t)], the sum extending over all the 
particles of the species a G {A,B}. Fig. [1] shows, for a low temperature T = 0.4, 
results for A- A, B-B, and A-B pairs. Intrachain static structure factors (i.e., chain form 
factors), S^ m (q), are also displayed. The latter quantities are computed by restricting 
the product p a (q, 0)p* (g, 0) = J2j k exp{iq • [r°(0) — r£(0)]} over pairs of monomers j, k 
belonging to a same chain. A sharp maximum is observed in Saa(<7) at q — 4.5, which 
corresponds to a typical distance of 1.4 between A-monomers. Results for Saa(<?) are 
qualitatively similar to those reported for the homopolymer case in a similar bead-spring 
model [29], [31]. A weak low-g structure is observed in the present case, which originates 
from the presence of 'holes' in the matrix of A-monomers. These holes are created by 
the inclusion of the B-monomers. The negative values of Sab (?) observed at small wave 
vectors are a signature of anticorrelation effects between A- and B-monomers at large 
distances, and indicate a moderate degree of demixing. This feature is illustrated in Fig. 
[21 which shows a typical configuration of the B-chains. The latter are not homogeneously 
distributed but form a sort of cluster structure. 

The partial static structure factor for B-B pairs, Sbb^), exhibits a rather different 
g-dependence (Fig. [T]). From a comparison with the form factor for B-chains, S^ m (q), 
it is clear that Sbb^) is largely dominated by intrachain contributions, as expected for 
high dilution of the B-chains in the matrix formed by the A-chains. The peak at q = 7.2 
corresponds to a typical distance of 0.87 between B-monomers. 

Data for the chain form factors in Fig. [1] are also compared with the Debye 
function [321 E3], S^ bye (q) = 2iVg- 4 (i^)- 4 {exp[-g 2 (i^) 2 ] + q 2 (R^) 2 - 1}, which is 
obtained by assuming a Gaussian distribution of monomer-monomer distances within 
the chain [321 ES]- As previously observed for the homopolymer case [3T], Gaussian 
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Figure 1. Partial static structure factors, <Saa(<7), Sbb{q) and Sab(q), at T = 
0.4. Also included are the chain form factors, Sf^ m (q), S B i B m (q), as well as the 
corresponding Debye functions, 5^ A ye (q), S BB hyc (q). Arrows indicate the wave vectors 
q = 2n/R" g , where and i?g are respectively the chain end-to-end distance and 
gyration radius of the species a. 




Figure 2. A typical configuration of the B-chains. 



statistics approximately work at low q but clearly break down for wave vectors probing 
distances smaller than the chain gyration radius . The magnitude of the deviations of 
simulation data from the Debye function is similar to observations for the homopolymer 
case [31J. Hence, chain statistics is not significantly affected by blending. 
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4. Main predictions of MCT 

In this section we summarize some of the main predictions of the Mode Coupling 
Theory (MCT) for the glass transition. Extensive reviews can be found, e.g., in Refs. 
[j~5| fi~6l [341 l35| [36] . In its ideal version, MCT predicts a sharp transition from an 
ergodic liquid to a non-ergodic glassy state for a given value of the relevant control 
parameter £ (in the following the temperature, T, though results exposed below are 
valid for any £). On approaching the transition point T = T c , MCT establishes a 
set of quantitative predictions for any correlator coupled to density fluctuations, $(£). 
An example is normalized density- density correlators of wave vector q, F aa (q,t) = 
(p a (q,t)pa(q,0)) / S aa (q), for the species a. At the critical temperature T = T c the 
long-time limit of <&(£) jumps from zero to a non-zero value, denoted as the critical non- 
ergodicity parameter, $ c . In the standard case (type-B transitions) the jump in $(£) is 
discontinuous, i.e., <3> c takes a finite positive value. 

For ergodic states close to the transition point $ usually exhibits a first decay to 
a plateau, whose time extension increases as the transition is approached. This plateau 
regime corresponds to the temporary trapping of each particle within the cage formed 
by its neighbouring ones, i.e., the well-known caging effect which is generally present 
in, e.g., supercooled liquids or jammed systems. At times longer than the so-called first 
MCT time scale t a , the correlator $ starts a second decay from the plateau to zero. This 
second decay is commonly known as the a-process and represents the full decorrelation 
of the system from its initial configuration, i.e. the structural relaxation. According to 
MCT, the initial part of the a-process (denoted as the von Schweidler regime) is given 
by a power law decay oc — t b , with < b < 1. A power-law series expansion extends the 
description of the a-decay to longer times: 



The prefactors h$ and h$ are state point-independent and are different for each 
correlator $. On the contrary, the von Schweidler exponent b is common to all 
correlators. The characteristic time scale of the a-relaxation, r a , is the second MCT 
time scale. It is also unique for all correlators, and diverges at the transition point 
as oc (T — T c )~ 7 (see below). The a-decay can often be described by an empirical 
Kohlrausch- Williams- Watt (KWW) function, oc exp[— (t/r)^*], with a ^-dependent 
stretching exponent < /3$ < 1. An interesting prediction of MCT [37J is that 
(3 q = b and r oc q~ l / h in the limit of large q, both for density-density [F(q,t)\, 
and self-correlators, F*(q,t) = exp{iq • [r a j(t) — r aj:7 -(0)]}\ /(NN a ). This result 

[231 EHl EHl EH EE1 ESI SQl IS] provides a consistency test for data analysis. 

Another prediction of MCT for state points close to the transition point is the 
power law dependence of the diffusivity and the relaxation time r*: 



The relaxation time r* of the correlator $ is defined as the time where <&(£) decays 
to some small value x, provided it is well below the plateau. The time-temperature 



= $c _ h ^t/ Ta f + hf{t/r a f b + 0(t 3b ). 



(1) 



r *j D cx (T — T c ) 



(2) 
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superposition principle of the MCT establishes that, for t much longer than the first 
time scale t a , $(t/r Q ) = $, where $ is a ^-dependent scaling function. According 
to this prediction, for a fixed x the ratio t*/t q is temperature-independent for any 
t* 3> t a , i.e, r* oc r a . In other words, r* will be ^-modulated but will follow the same 
power law behaviour in T as the a-relaxation time r a [even if V* r a ). Note that, in 
the MCT terminology, r a is a single time scale, though its value can be approximately 
probed by evaluating dynamic correlators $ for which r* ~ r a . This is the case of, 
e.g., the density-density correlator F(q,t) for wave vector q at the maximum of the 
static structure factor S(q), since the former probes decorrelation over typical distances 
between nearest-neighbour particles. The relaxation time of F(q,t) is indeed often 
denoted as the 'a-relaxation time', though in the context of MCT the latter strictly 
corresponds to r a . 

The exponent 7 in Eq. (j2J) is given by the relation: 

with < a < 0.395. Hence 7 > 1.766. The critical exponents a, b, and 7 are univocally 
related with the so-called exponent parameter A through: 

A = r2 ( 1+6 ) = r2 ( X - a ) (A) 
r(l + 26) r(l-2a)' 1 1 

where V is the Gamma function. The exponent parameter A is univocally determined 
by the static correlations (i.e., by the total and partial static structure factors) at the 
transition point T = T c . For type-B transitions it takes values 1/2 < A < 1. 

When numerical solutions of the MCT equations are not available the non- 
ergodicity parameters, prefactors and exponents in Eqs. ( jT1T2"1l3~f4l) — which are system- 
dependent quantities controlled by static correlations — are empirically obtained as 
fit parameters from simulation or experimental data. Consistency of the data analysis 
requires that the so-obtained set of exponents fulfill both Eqs. (J3j) and (BJ. 



5. Dynamics of the slow component in the blend 

Fig. [3] shows results for the mean squared displacement averaged over all the monomers, 
(Ar^(t)), both for A- and B-chains. The introduction of monomer size disparity, 
oaa/obb = 1.6, induces a large time scale separation between the two components, 
for low concentration of the B-chains, by decreasing temperature. Now we analyze 
relaxation features for the slow A-component. Results for the fast B-component are 
analyzed in the next section. 

Figs. EJ EJ [6j [7] and [8] show, for the A-chains, a consistent test - - i.e., with 
a common set of exponents — of MCT predictions for several dynamic correlators, 
diffusivities, and relaxation times. Fig. H] shows, for several wave vectors, results at 
T = 0.4 for the intrachain coherent correlator F^ m (q,t). The latter is computed as 
( X^-.fc^P^q ' [^j'W - r a,fc(0)]}^) /[NN a S^ m (q)], for any species a. In this equation 
the sum only includes j, k pairs belonging to a same a-chain. Fig. [5] shows normalized 
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Figure 4. Symbols: For different wave vectors, intrachain coherent correlator for A-A 
pairs, F^ in (q, t), at T = 0.4 . Lines are fits to Eq. (P) with an exponent b = 0.30. 



correlators of the Rouse modes 4>pp(t) at T = 0.5. The latter are defined as 
= (X£(*) • X£(0))/([X^(0)] 2 ), where the Rouse normal modes [321 [33] of index 
p — 0,1, ...N — 1 are given by X p (i) = iV _1 X^=i r j(^) cos[jpn/N]. Fig. [6^ displays, at 
T = 0.45, angular correlators C^ A (t) for the bond vector, b(t), between consecutive 
monomers. Such correlators are defined as C^ A (t) = P n [cos6(t)], where P n is the 
Legendre polynomial of order n, and cos6(t) = (b(t) -b(0))/ (b 2 (0)). Angular correlators 
Cn^it) for the chain end-to-end vector, e(t), are defined in analogous way, with 
cos 6(t) = (e(t) ■ e(0))/(e 2 (0)). Data for cL° )A (t) at T = 0.5 are given in Fig. [6b. 

Lines in Figs. HI [3, and [6] are fits of the a-decay of the mentioned correlators to a 
power-law series expansion as Eq. (Tj[|) with a common von Schweidler exponent b = 0.30. 
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Figure 5. Symbols: For different values of p, correlators of the Rouse pth-modes, 
4>pp(t), of the A-chains. Lines are fits to Eq. fl} with an exponent b = 0.30. The 
temperature is T = 0.5. 



(a) 




Figure 6. Panel (a): Symbols correspond, for different values of n, to angular 
correlators C^ A (t) of the bond vector of the A-chains. Lines are fits to Eq. (fT]) 
with an exponent b = 0.30. The temperature is T = 0.45. Panel (b): As panel (a) for 
angular correlators C„ (t) of the end-to-end vector. The temperature is T = 0.5. 

Only terms up to second order (t 2b ) are included in the fit procedure (in the following, 
references to this equation will be understood as limited to second order). It must be 
stressed that the validity of Eq. (CQ) for the early-middle a-decay must not be assessed by 
the length of the vertical interval of that it is able to cover. Indeed, if the relaxation 
time of the analyzed correlator is much longer than the ct-time r a , the vertical interval 
described by (CQ) will be rather small, as we will discuss below. The prefactors h$ and 
h<$> m CQ) , which yield the amplitude of the decay, are generally in anti-phase with $ c 
[Ml 1361 1381 W2\ and are small for large values of the latter. Hence, for correlators with high 
plateaux Eq. ([!]) will only describe a small vertical interval of the decay. On the contrary, 
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Figure 7. Panel (a): (/-dependence of the stretching exponent of different correlators 
for the A-monomers (see text for notations) at temperature T = 0.4. The dashed 
line indicates the large-g limit (3(q) = b = 0.30. Panel (b): As panel (a) for the 
corresponding KWW times (see text). The dashed line corresponds to the power law 
cx q-^ b , with b = 0.30. 



validity of ([I]) is given by the extension of the time window (i.e., horizontal interval) 
that it is able to describe. In the present case a good description of the simulation 
data is obtained over three time decades for the lowest investigated temperature, a time 
window of validity which is typically achieved in simulations. It must be noted that 
such a time window corresponds to a specific dynamic regime, the early-middle stage of 
the structural a-relaxation. However, relaxation of a given correlator to a small value 
(e.g., x = 0.2) can occur at a very different time scale r*. This is the case of, e.g., 
low-index correlators of Rouse modes or chain end-to-end vectors. The latter show a 
decay much slower than density-density correlators -Faa(<?, t) at the maximum of Saa(<z) 
(q = 4.5), which properly probe the time scale r Q of the structural a-relaxation for the 
A-chains. For Fj^(q,t) we find To. 2 = 1.7 x 10 4 at T = 0.5, while relaxation times at 
the same temperature for low indexes of &p P (t) and Cn (f) are clearly much longer (see 
Figs. [5] and [6b)- Having said so, MCT establishes than asymptotic expansions as (pQ) 
will be observed for any dynamic correlator in the specific early-middle time window of 
the a-relaxation, the proccess here investigated. 

Fig. [7] shows a test of the MCT predictions f3(q — > 00) = b, and r(q — > 00) oc q~ x l h . 
The stretching exponents (3, /3 cha , and /3 s , correspond respectively to the density- density, 
intrachain coherent and self-correlators of the A-monomers at T = 0.4, and are obtained 
as fits of the decay from the plateau to a KWW function. The corresponding KWW 
times are respectively denoted as r, r cha , and r s . The mentioned large-g predictions 
for stretching exponents and KWW times are fulfilled with b = 0.30, i.e., with the 
same value of the von Schweidler exponent used in the fits of the dynamic correlators 
presented in Figs. HI El and El 

The exponent b = 0.30 provides through Eqs. (j3J) and (j4j) the values A = 0.90, 
a = 0.21, and 7 = 4.0. Now we test the validity of Eq. (j2J) with this latter value of 
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Figure 8. Symbols: Inverse diffusivity (circles) and relaxation times of different 
correlators for the A-monomers (see text for notations). The wave vector for To. 2, 
Tg^ ain , and Tq 2 is q = 4.6. The dashed lines are (from top to bottom) fits of £^cm> r^f, 
r s 2 , and r b | to the MCT power law oc (T — T c )-T, with T c = 0.37 and 7 = 4.0. 



7. Fig. M shows the temperature dependence of the relaxation times r* of several 
dynamic correlators $ for the A-monomers. As mentioned above, these times are 
defined as those where the corresponding correlator decays to a value x. Notations T0.2, 
TQ 2 am , Tq 2 , To!, Tq|, Tq 1 !, and r^f correspond, respectively, to the correlators F\j^(q,t), 
F^ in (q,t), Fl(q,t), Cf )A (t), Cf )A (t), ^ 8 (t), and ^(t). The wave vector for the first 
three correlators is q = 4.6, an intermediate value between the main maxima of Saa{q) 
and Sfj^ n (q) (see Fig. [Q. The selected values of x are well below the plateau height of 
the correlator (see Figs. HI [U and [6]). Also displayed is the inverse diffusivity, -D^M' °f 
the center-of-mass of the A-chains, which is defined as the long time limit of the ratio 
6t/([Ar^ M (t)] 2 ), where ([Ar^ M (t)] 2 ) is the corresponding mean squared displacement. 

Dashed lines in Fig. [8] represent fits to the power law (j2J) by forcing a common 
critical temperature T c for all the data sets, with a fixed exponent 7 = 4.0, i.e., the 
value independently determined from the analysis previously presented in Figs. HI El E] 
and [71 A value T c = 0.37 provides the best global fit with the mentioned constraint. 
Interestingly, this value is much lower than the one obtained for the homopolymer state 
[14"] at the same investigated packing fraction = 0.53, T c = 0.52. The latter value is 
obviously identical for the limits x-q = and x-q = 1, since the energy scale e of the 
model is the same for all the pair interactions, which only differ by the length scale 
a a/ 3 (see Section 2). Hence, blending at a fixed packing fraction stabilizes the ergodic 
phase as compared to the homopolymer state, in analogy with the behaviour observed 
for colloidal binary mixtures of similar size disparity [121 H3j HI] . 

A good description of all the data sets is obtained over more than two decades in 
relaxation time and inverse diffusivity As expected, due to the asymptotic character 
of Eq. (j2]), deviations from power-law behaviour occur at high temperature. Such 
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deviations are also present below some ill-defined temperature very close to T c . This 
feature is often observed if one investigates dynamics at sufficiently low temperatures 
[4"5[ SHI HZ! an d is usually related with the presence of activated hopping events, which 
are not accounted for within the ideal version of the MCT. 

The analysis of simulation data of the slow A-component that has been presented 
in Figs. HI and [8] consists of a series of independent tests of several predictions 
of MCT with a common set of values of the critical exponents. Therefore it provides 
a robust determination of such values, and in particular of the exponent parameter 
A = 0.90 from which the rest of the exponents are derived through Eqs. ([3]) and (J4j). This 
value of A is unusually large, as compared to those typical of one-component systems, 
as monodisperse hard spheres [36J (A = 0.74), simplified models of orthoterphenyl [39] 
(A = 0.76), silica [48] (A = 0.71), water [49] (A = 0.78), or bead-spring homopolymers 
[50] (A = 0.72). In the following subsection we discuss the consequences of the large 
value of A here obtained. 

6. Dynamics of the fast component in the blend 

Now we analyze the dynamics of the fast B-component. Fig. [9h shows simulation 
results for the intrachain coherent correlator, F^ m (q,t), at temperature T = 0.4. As 
previously reported in Ref. [H] for the total density- density correlator F BB (q,t), a 
concave-to-convex crossover is observed by varying the wave vector. For intermediate 
values of the latter, a purely logarithmic decay occurs over more than three time decades. 
Following a procedure analogous to that of Ref. [T4] , we have analyzed the decay of 
Fgg am (g,t) in terms of a logarithmic expansion, 

^BB ail %, t) = f c q - H q \n{t/r a ) + \n\t/r a ) + 0[\n 3 (t/r a )], (5) 

with To- ~ t a , instead of the von Schweidler series (CD) used for the A- monomers (we will 
discuss this point below). Within the framework of MCT, logarithmic expansions of 
dynamic correlators are associated to the presence of a nearby higher-order transition 
[T7] [TBI [T9~] . The latter is characterized by a value of the exponent parameter A = 1, 
though analogous predictions are expected for sufficiently large values A — > l - as the 
one here obtained, A = 0.90 (see Section [5]). Eq. ([5]) provides a good description of 
the decay of correlators displayed in Fig. [9^. Analogous fits are shown in Fig. [9b for 
the orientational correlators C„ B (t) of the bond vector, evaluated for different values 
of n. In this case the validity of the logarithmic expansion is observed, at the same 
temperature, over a shorter time interval. Fig. [10] shows the values of the coefficients 
fq, H q , and H^p obtained from the corresponding fits of F B i B am (q,t) at two different 
temperatures (T = 0.4 and T = 0.5). The term f£ is the critical non-ergodicity 
parameter, which is associated to the transition point. Therefore its values at different 
wave vectors must not depend on the state point at which they are obtained as fit 
parameters. This is confirmed by the numerical values displayed in Fig. [TUh . According 
to MCT, the prefactor H q is factorized as the product of two terms. One of them only 
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Figure 9. Panel (a): Symbols correspond, for different wave vectors, to intrachain 
coherent correlators for B-B pairs, F-^ m (q,t), at T = 0.4. Lines are fits to the 
logarithmic expansion ([5|). Panel (b): As panel (a) for the angular correlators C^ B (t) 
of the bond vector of the B-chains. 

depends on the state point and the other one on the wave vector [19J. Therefore the 
values of H q evaluated at different state points must obey scaling behaviour. This feature 
is also confirmed by data in Fig. [TUb . Also in agreement with MCT expectations [19J, 

(2) 

the obtained values of the second prefactor H q are smaller than H q and uncompatible 
with scaling behaviour (see Fig. [TUb). 

Results presented in Figs. [9j and [10J support the similar analysis performed in Ref. 
[H] for density-density correlations of the B-monomers, Fbb(?, t). It must be stressed 
that the choice of Eq. (J5J) for describing relaxation of correlators for B-monomers is 
not, in principle, in contradiction with the description of the same correlators for the 
A-monomers in terms of the power-law series (CD). Both equations are series expansions 
whose convergence depends on the analyzed region of the control parameter space. For 
the case of higher-order transitions (A = 1), or more generally for transitions with 
A — > l - , there are g-dependent paths in the control parameter space where the series 
(JSJ) is rapidly convergent. In particular for each wave vector there are optimal paths 

(2) 

where Hq = 0. Along these paths the corresponding correlator will exhibit a purely 
logarithmic decay [19J. Moreover, by properly tuning the control parameters or the 

(2) 

wave vector, it is possible to change the signus of H q and, as a consequence, inducing 
a concave-to-convex crossover in the shape of the decay |19j . as observed in Fig. [9j 
Since from the analysis of dynamic correlators for the A-monomers we have determined 
a value A = 0.90, it might be expected that such correlators will exhibit such features 
at some state point. Indeed, they are observed at higher temperatures, as shown in Fig. 
[U]for Fp^ m (q,t) at T = 1.0. The decay exhibits a clear concave-to-convex crossover by 
tuning the wave vector. Logarithmic relaxation covers two time decades for q rs 5.2. 

The fact that features associated to nearby higher-order MCT transitions are 
observed for the A- and the B-component at very different temperatures must be 
commented. As mentioned above, the optimal paths in the control parameter space 
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Figure 10. Symbols in panels (a), (b), and (c) correspond respectively to the values 

(2) 

of the critical non-ergodicity parameter, /£, and the prefactors H q and H q ' in Eq. 
©, for the intrachain coherent correlator for B-B pairs, F^ m (q, t). Temperatures are 
T = 0.4 (circles) and T = 0.5 (squares). 




Figure 11. Symbols: For different wave vectors, intrachain coherent correlator for 
A-A pairs, F^ in (q, t), at T = 1.0. The straight line indicates logarithmic behaviour 
over two time decades. 
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Figure 12. Comparison, at T = 0.4, between the intrachain coherent (Fj^ m (q,t), 
symbols) and self- (F^(q,t), dashed lines) correlators for the B-monomers. The wave 
vectors for the latters are, from top to bottom, the same ones as for the formers (see 
legend). 



for the observation of logarithmic relaxation are different for each correlator [19] . The 
location of these paths is controlled by static correlations [19], which in the present 
case are very different for the A- and the B-monomers (see Fig. [1]). This difference 
might explain why anomalous relaxation for both types of monomers is observed at 
very different temperatures. Unless one moves close to the optimal path in the control 
parameter space — which usually involves a simultaneous variation of several control 
parameters [I7J [19] — logarithmic relaxation vanishes by decreasing temperature, and a 
standard two-step decay is recovered (see Ref. [21] for an illustrative example). This 
seems to be the case of correlators for the A-monomers, which at low temperature are 
well described by the von Schweidler power-law series (0Q). A similar result has also been 
observed for mixtures of large and small non-bonded particles [23] . Still, a satisfactory 
answer to this point can only be obtained by solving the MCT equations for this system. 

Fig{12] displays for the B-chains at T = 0.4, results for the density self-correlators, 
F#(q,t). A reliable fit (over more than one time decade) of the corresponding decays to 
Eq. (0Q) — with any common exponent b for all wave vectors — or §5§ was not possible. 
Similar tests were also unsuccesful for correlators probing reorientations of chain end- 
to-end vectors, Cn^ B (t), and relaxation of Rouse modes, of the B-chains. The 
reason for the apparent failure for the B-chains, or at least limited range of validity, of the 
former equations for these correlators remains to be understood. It might be that this 
feature is connected to a non- universal character of the asymptotic expansions (pQ), ©for 
binary mixtures with very different time scales for their respective density fluctuations, 
and that despite this non-universality, MCT can still reproduce the behaviour of the 
mentioned correlators for the B-chains. It might either be related to the presence of 
hopping events intervening in self-motions of the B-monomers (see below) . As we argue 
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in the following, there are results in the literature that support these possibilities. 

Numerical solutions of the MCT equations have been recently presented for sodium 
silicate melts and compared with simulation results [51]. In these systems the fast 
sodium atoms and the slow silica matrix exhibit a strong time scale separation similar 
to that observed here for the A- and B-chains [511 E2j. Though an analysis of density- 
density correlators for the different atomic species, as well as of self-correlators for silicon 
and oxygen, have provided a consistent test of MCT predictions with a common set of 
dynamic exponents [52], self-correlators F s (q,t) of the sodium atoms do not show |52j, 
as in the present reliable time interval for apparent validity of Eqs. ([I]) or ([5]). 

Still, the corresponding numerical solutions of MCT equations reported in Ref. [51] 
do reproduce the qualitative behaviour of F s (q,t) for the sodium atoms. In particular, 
MCT gives account for the unusual time scale separation between self- and collective 
density correlators which is observed for the sodium atoms. This feature is assigned 
[52| 1531 EU [551 EH] to preferential diffusion along a long-living structure of channels 
induced by the much slower relaxation of the silica matrix, which leads, for the alkali 
ions, to a fast decay of self-correlations as compared to collective correlations. 

Fig. [12] shows a comparison between self- and intrachain coherent correlators for 
the B-chains at T = 0.4. Both correlators only converge to each other in the limit of 
large-g. Since, due to the low concentration of the B-component, intrachain coherent 
correlators for the B-chains exhibit only small differences (not shown) with density- 
density correlators for all the B-B pairs, the large time scale separation between F%(q, t) 
and F^(q,t) presented in Fig. [12] is a feature analogous to that above commented for 
alkali ions in silica matrices [52]. Indeed, following a procedure similar to that presented 
in Ref. [53] for the sodium atoms, we have determined a similar structure of channels 
for preferential motion of the B-chains. We have divided the simulation box in cubic 
subcells of size ~ a BB and computed, for a trajectory of the system, the number of 
times each subcell is visited by a B-monomer. Fig. [13] displays, at T = 0.5, the N&N 
(a number equal to that of B-monomers) most visited subcells for a simulation time 
t = 5 x 10 4 . The latter is much longer than the time for structural relaxation of the 
B-monomers at that temperature As shown in Fig. [131 the mentioned subcells are not 
randomly distributed but form connected clusters, in analogy with results reported in 
[53] for sodium in silica matrices. 

In Fig. [13] we also display the initial and final configuration of the B-monomers 
for the mentioned simulation interval t = 5x 10 4 used for the computation of the most 
visited subcells. As expected (t is much longer than the structural relaxation time for B- 
monomers) both configurations are fully decorrelated. Therefore the mentioned channel 
structure is not a trivial consequence of the static correlations for the B-B pairs, which 
also form a cluster structure (Fig. [2]). It is instead induced by the time scale separation 
of the dynamic correlations, which are much slower for the confining matrix formed by 
the A-chains. The channel structure will only vanish when any region of the simulation 
cell will be visited by the B-monomers with the same probability. This can only occur 
at much longer times probing full structural relaxation of the A-component. A detailed 




Figure 13. Left side: Cubic boxes represent the N-gN subcells (of size rj 0bb) in 
the simulation box which have been visited more times by the B-monomers during a 
simulation time of t = 5 x 10 4 , for T = 0.5. Right side: Initial (dark spheres) and final 
(light spheres) configuration of the B-monomers for the latter simulation interval. The 
same orientation of the simulation cell is used in both figures. 



static and dynamic characterization of this channel structure is beyond the scope of this 
article and will be presented elsewhere. 

Finally, it must be mentioned that the observed decoupling between intrachain 
collective and self-correlators is exhibited only by the B-chains in the blend. For 
the A-chains in the blend, as well as for the homopolymers, we have observed only 
small differences for the latter correlators. Decoupling between self- and collective 
intrachain dynamics is indeed a rather unusual feature, at odds with expectations from 
the standard Rouse model [32j [33]. This observation in the simple bead-spring blend 
here investigated is supported by recent neutron scattering experiments on PEO /PMMA 
at low concentration of PEO [57]. Whether numerical solutions of MCT equations are 
also able, in analogy with the case of alkali ions in silica, to give account for this feature 
is an opened question. 

Fig. E] shows the temperature dependence of the relaxation times for several 
dynamic correlators probing relaxation of the B-component. Notations r 02 , ^o^ 111 , 
Tq 03 , Tq\, Tq 2, and correspond, respectively, to the correlators F BB (q,t), F B i Q aATL (q,t), 
F B (q,t), C[ (t), C[ (t), and <f>ii(t). The wave vector for the first three correlators 
is q = 4.6. Also included is the inverse diffusivity of the center-of-mass for the B- 
chains. The set of data shown in Fig. [TH exhibits a behaviour rather different from 
similar quantities for the A-chains displayed in Fig. [HJ Only relaxation times for 
collective density correlations, F BB (q, t) and FgB am (g,t), show qualitative agreement 
with the MCT power-law oc (T — 0.37)~ 40 derived from data of the A-monomers. As 
expected, deviations occur at temperatures very close to T c . The rest of the quantities 




Figure 14. Symbols: Inverse diffusivity (circles) and relaxation times of different 
correlators for the B-monomers (see text for notations). The wave vector for To. 2, 
T o h i3 m > and Tq 03 is q — 4.6. The thick and thin solid lines are fits of, respectively, r .2 
and r^ n to a MCT power law tx (T - T c )-^, with T c = 0.37 and 7 = 4.0. The arrow 
indicates the inverse value of T c . The dashed lines are (from top to bottom) fits of 
-^CM' t 02j t o2i an d Tq 03 to Arrhenius behavior, oc exp(E/T). The activation energies 
are respectively E — 3.4, 4.1, 3.8, and 5.9. 

displayed in Fig. O are uncompatible with power-law behaviour. They show instead an 
apparent Arrhenius dependence, oc exp(E/T), from moderate to the lowest investigated 
temperature. This feature is demostrated in Fig. [141 by the linear behaviour observed by 
representing data in logarithmic scale vs. (linear) 1/T. The obtained activation energies 
E vary between 3.4 for the center-of-mass diffusivity and 5.9 for the relaxation time of 
density self-correlations. The observed Arrhenius behaviour suggest that strong hopping 
events intervene in the structural relaxation of the B-chains, similarly to observations 
for alkalin ions in silica [58]. These events seem to affect more strongly to self- than to 
collective density correlations, for which a power law behaviour can be observed over 
two decades in relaxation time for temperatures above T c . It remains to be understood 
whether such hopping events — which are not included in the ideal version of MCT - 
are related to the mentioned reduction of the range of validity of Eqs. (TO or (151) for the 
corresponding correlators. It is worth mentioning that the latter possibility might be 
the case for sodium atoms in silica. Numerical solutions of the MCT equations reported 
in Ref. [5T] . though reproducing the observed qualitative behaviour, understimate the 
strength of the decay exhibited in simulations for self-correlators F s (q,t) of sodium. 
Hence, the presence of hopping events presumably accelarates relaxation as compared 
to theoretical predictions. 

Finally, it is worthy of remark that the observed Arrhenius- like temperature 
dependence for Tq % and Tq 03 is consistent with experimental observations, for the fast 
component, in real polymer blends with large dynamic asymmetry by, respectively, 
dielectric spectroscopy [3, EJ [TO] and neutron scattering [15] . which probe relaxation 
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times of similar dynamic correlators. Arrhenius behaviour for self-dynamics is also 
observed for the case of alkali ions in silica [58] or for water reorientation in polymer 
matrices [59]. This common Arrhenius-like behaviour in very different systems suggests 
a universal feature for low concentrations of fast molecules in slow host media with 
interconnected voids. 

7. Conclusions 

We have presented a computational investigation on the structural relaxation of a 
simple bead-spring model for polymer blends with large dynamic asymmetry. We have 
computed a large set of dynamic correlators probing relaxation of density fluctuations, 
Rouse modes, and reorientation of bond and chain end-to-end vectors. Results have 
been discussed within the framework of the Mode Coupling Theory (MCT) for the glass 
transition. A robust test of MCT predictions has been achieved through a description 
of the different analyzed correlators with a common set of dynamic exponents, though 
for some correlators probing dynamics of the fast component MCT asymptotic laws are 
apparently not observed. The observation of Arrhenius-like behaviour suggests that this 
breakdown might be associated to strong hopping events intervening in relaxation of the 
fast component. 

An unusually large value of the exponent parameter A has been obtained, close to 
the upper limit (A = 1) characteristic of higher-order MCT transitions. According to 
MCT predictions, the anomalous relaxation features observed in the present system, 
as logarithmic decays or concave-to-convex crossovers in density correlators, might be 
associated to that underlying higher-order scenario. An investigation of the case of 
extreme dilution, where each individual chain of the fast component is sorrounded 
only by chains of the slow component (and where an asymptotic dynamic limit is 
expected [6]), would be computationaly expensive. Still, we expect that a qualitatively 
similar scenario of anomalous relaxation will be observed. Since chain connectivity will 
always guarantee the presence of neighbouring monomers of the same species for a given 
monomer of the fast component, coexistence of bulk-like caging and confinement for the 
fast component would be present even at extreme dilution, inducing the higher-order 
scenario. On the other hand, a progressive increase of the concentration of the fast 
component will reduce the time scale separation (i.e., the dynamic asymmetry) between 
the two components, and confinement effects will finally vanish. In that situation a 
standard MCT relaxation scenario (as observed for the homopolymer case) will be 
recovered. The large collection of results here presented might motivate theoretical 
work on structural relaxation in polymer blends with large dynamic asymmetry, and in 
particular, numerical solutions of the MCT equations to confirm the suggested higher- 
order scenario. 
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